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Abstract 

Unparticles (JA) interact weakly with particles. The direct signature of unparticles will be in 
the form of missing energy. We study constraints on unparticle interactions using totally invisible 
decay modes of vector quarkonia V and neutrinos. The constraints on the unparticle interaction 
scale are very sensitive to the dimension dy of the unparticles. From invisible Z and V decays, 
we find that with djj close to 1 for vector lA, the unparticle scale Ky can be more than 10^ TeV, 
and for dy around 2, the scale can be lower than one TeV. From invisible neutrino decays, we find 
that if dy is close to 3/2, the scale can be more than the Planck mass, but with dy around 2 the 
scale can be as low as a few hundred GeV. We also study the possibility of using V{Z) ^ + U 
to constrain unparticle interactions, and find that present data give weak constraints. 
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Introduction 

Recently Georgi proposed an interesting idea to describe possible scale invariant effects 
at low ene.g.es by a„ operate. 0„, te.„.ed unparticle Q. Based on a specific scale invariant 

theory with a non-trivial infrared fixed point by Banks and Zaks [2], it was argued that 
operators Obz made of BZ fields may interact with operators Osm made of Standard Model 
(SM) fields at some high energy scale by exchanging particles with large masses, Mk, and 
induce interactions of the form 

-OsmObz, (1) 



where cIbz and dsM are the dimensions of the operators Obz and Osm- 

At another scale An the BZ sector induces dimensional transmutation, and below that 
scale the BZ operator Obz matches on to unparticle operator Ou with dimension djj. The 
unparticle interaction with SM particles at low energy then has the form 

/A \ dsM+dBZ-4: 

CuXAtr''"-'''OsMOu, A=(^j . (2) 

The unparticle may have different Lorentz structures such as a scalar Oy, a vector Oy, 
a spinor O^, and etc.. The specific form of SM particle and unparticle interactions are 



not known and are usually parameterized in terms of operators. In Ref. [17|] a class of 
operators involving SM particles and unparticles are listed. Using these operators one can 
study unparticle phenomenology in a systematic way. 

One of the main phenomenological goals of unparticle physics study is to find out at 
what energy scale unparticle effects may show up [1, 3-27]. The most direct signatures 
of unparticles will be in the form of missing energy in invisible decays of particles with 
an unparticle U in the final state. In this paper we study constraints on the unparticle 
interactions using Z invisible, V invisible, and v —* invisible decays. We will also 
study the possibility of using V{Z) —>■'-/ + to constrain unparticle interactions. 
Constraints from invisible decay of Z boson 

The process Z ^ U contributes to invisible decay of Z. For this process the following 
operators, with SM fields and derivatives contribute less than or equal to 4 dimensions, will 



contribute [Jj 

Ko^l-'-B,^d^Ol, , KoAl'^B^^d^Ol, , X',,Alr'nH^D,H)0^ . (3) 
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Here the vector unparticle operator Ojj is hermitian and transverse with, 9^0^ = 0. 
The matrix elements for Z resulting from the above operators are given by 

M{Z ^ \',o) = Ko^u"" sinOwikz ■ koez ■ eo - kz ■ eoko ■ ez), 
M{Z ^ U^, ~X',o) = iKo^u"'' sinewe^uaf^k^z^^zkyo, 

M{Z A',,) = -2/rn(A;,,)A^-'^"— ■ eo, (4) 

Sm(^/(7vy ) I 

where v = 246 GeV is the vacuum expectation value of the Higgs doublet. Since kz = ko 
for Z to transition, the second term in the above does not contribute. Here we have 
used the notation ez and ep to describe the polarizations of Z and U^. 

For a decay of a particle into an unparticle and other particles, the differential decay rate 
is given by Ij 

dTiP^l()=^p^d^P), (5) 
2mp 

where (i$(P) is the phase space factor for the decay. It is given by 

rf$ = Ji2nmp-Y.p,)Y[d^p,)^^. (6) 

For a particle the phase factor d^{pj) is equal to 27r^(p°)(5(pj — m|), and for an unparticle it is 
given by AdJ{py{p^){py^'^ with A^^ = (167r5/V(27r)2'^")(r(di, + l/2)/r{du - l)r{2du)). 
The decay width of a particle decay into an unparticle is given by 

r(P-w) = ^^^J^p)'"-'- (7) 

We note that in the limiting case of rf^^ = 1 the decay width becomes zero since A^i^ 
has a factor l/r{d — 1) which goes to zero when du 1. Physically this is because that 
in this case jll, \imd^^i+ Ad^^Oip^) / p^^'^~'^^'i = 27i6{p'^), the unparticle behaves as a massless 
particle. When 7^ 0, the decay rate T{P — > U) is zero. 

For Z U, we have 

r(z ^ Aio) = ^ sin' e„(y„fUi.\ a,^ , 



In general four parameters are needed to describe unparticle interactions with SM par- 
ticles: Cu, A, Ku and du as shown in eqs.([T]) and ([2]). If dsu + dsz > 4, the parameters 



A is less than 1. The parameter Cu contains information about the original heavy particle 
mediating interaction of the SM and scale invariant sectors, and also information about 
the transmutation. One may normalize the parameter Cy into the definition of A for one 
operator, but in general will not be able to do so for more than one operator. The values 
for Cjj depend on the detailed dynamics. If one is only concerned with the scale where dif- 
ferent transitions have happened, one usually sets Cu to be one and use the two parameters 
A and h^u to describe the situation. In our numerical discussions, we will also follow this 
subscription. 

Precise experimental data have been obtained on Z decay widths 29|| with the invisible 
width to be: T{Z invisible) = 499.0 ± 1.5 MeV. This is to be compared with the width 
of 501.65 ±0.11 MeV from SM prediction for Z decy into neutrinos. New contribution to 
invisible Z decay is therefore constrained severely, basically need to be within the range 
of experimental error bar of order one MeV. In Fig. [1] we show constraints on unparticle 
interactions allowing the unparticle contribution to invisible Z decay to saturate 2a error of 
experimental data of 3 MeV. Numerically the bound on Im{X'i^f^) is about 5.5 times stronger 
than X'fjQ for given du and An. In Fig. [T]we only show constraint on X[q. Since for a given 
du, the scale Ak also depends on the parameter A[,q, one can view the constraints on X'j^q 
for a given An or on the scale Ak for a given X[q. In any case, from Fig{T]it is clear that the 
constraints are very sensitive to the dimension parameter du. If du is close to 1, for example 
with X'^Q = 1 and du to be 1.3, Au needs to be larger than 10^ TeV, but Au can be as low 
as one TeV for du = 2. If by some means the scale Au is known, for example Au = 1 TeV, 
we have the upper bounds X'j^q = 0.049 and 0.10 for the cases (ie^=1.3 and 1.5, respectively. 
Constraints from invisible decay of quarkonia 

For a vector quarkonium V decays into an unparticle U the following operators will 
contribute. 

X'qqAI'^'QlI.QlOI:, X:^uAI^'-Ur^,UrO^, X'^M^'-Dni.DnOj^. (9) 

The matrix elements for vector quarkonia and unparticle transition resulting from the 
above interactions can be written as the following 

MiV^W^A') = lyK-'^'iOMV) ■ (10) 

where for q being an up type quark. A' can be X'qq and X'^^ with Qg = 2/3, and for q being 
a down type quark. A' can be X'qq and A^^, with Qg = —1/3. 
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FIG. 1: Bounds on parameter space of Ku vs. X'^q allowing unparticle decay mode to saturate the 
difference of 3 MeV for invisible decay width of Z between SM prediction and experimental data. 
The lines from top to bottom are for du = 1.3, 1.5, 2.0. 



We obtain 



Br{V ^W,X') 



du-i 



Br{V fi+fi~) 327ra2Q2 yj^^ 



(11) 



The operator AJ,qA^ '^^ B^yd^O" also contributes to this process. We have 
Br{V 



2a 



A2 



(12) 



Br{y fi+fi-) 

In Fig. [2] we show constraints on unparticle interactions using experimental data [28| : 
Br{T — > invisible) < 2.5 x 10"'^. For this case Q = D = b and Qg = —1/3. In obtaining 
the constraints, we have neglected small contributions from T ^ z/z/ to T invisible decay 
width. In this case, for given du and A^, the constraints on Agg and X'j^j^ are the same, 
while the constraint on X[q is 5.6 times weaker than the bounds for Agg In Fig. [2], we 
show constraints on Agg and A^. Again we see that the bounds are very sensitive to 



the dimension du- If du is close to 1, with Agg = 1 and dt 



than 2 X 10^ TeV, but Ku can be as low as 400 GeV for di 



1.3, Ku needs to be larger 
2. If the scale Ku is set to 



be 1 TeV, we find the upper bounds X'qq^dd = 0.025 and 0.082 for the cases du=l.?) and 
1.5 respectively. Note that the constraint on A'f,^ obtained from T — W is weaker than that 
obtained from Z ^ U hy a, factor of 1.49(m^/mx)'^"~^, which is in the range 2.9 ~ 14.4 
when du is in the range of 1.3 ~ 2.0. 
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FIG. 2: Bounds on parameter space of Au vs. X'qq allowing unparticle decay mode to saturate 
the experimental data for invisible T decay. The lines from top to bottom are for du = 1.3, 1.5, 2.0. 

One can easily work out the case for invisible decay of J/ip by taking Q = c, U = c and 
Qq = 2/3. BES has accumulated more than a million J/ip, it would be interesting to see if 
these data when analyzed for invisible decay, a better constraint could be obtained. 
Constraints from invisible decay of neutrinos 

For an active neutrino decays into an unparticle, the following operator will contribute 

KAH'-^'^LlHOI, . (13) 

After the Higgs develops its vev, one obtains a transition matrix element between a 
neutrino and an unparticle 

Miu W) = XAT-'-ul^OI,. (14) 



This leads to 



.i|A,|2^ . (15) 

Using constraints from solor neutrino data [33i] T/m> 10~^s/eV on neutrino lifetime and 
mass ratio, one can obtain information about the unparticle interactions. Since the term 
generating this invisible neutrino decay is related to the Yukawa coupling, it is natural to 
have Xsv/\/2 to be of order the neutrino mass itself (one can easily converts into different 
normalization). In this case one would obtain 

ml ( < 1.3 X lO-^W. (16) 




FIG. 3: Allowed parameter space for Ku vs. du- 
Applying this formula to the solar neutrino with the constraint on its relevant mass to 



be larger than J Am'j.^i^j., we obtain the bound on Au as a function of du in Fig. [3] with 



the central value of Am^^,^^ = 8.0 x 10"^ eV^ |29|. It can be easily seen that the bound 
on An depends on du very sensitively. If one sets An to be less than the Planck scale 
mp = 1.22 X 10^^ GeV, the dimension du must be bigger than 1.6. 
Constraints from radiative — > 7 + invisible 

Now we study the possibility of using V{Z) ^ 7 + W to constrain the unparticle interac- 
tions. Using the general formula in eq.(j5]), we obtain the differential rate for V ^ 'j + U for 
a given matrix element M, 

dE, - 2^p^'''^^^^ ■ ^^^^ 

For Z(y) lA process, the unparticle U must be a vector type. For Z{y) ^ •y + U, the 
unparticle can be a scalar or a vector. Using this process, constraints on scalar unparticle 
interactions can also be obtained. We find that the following operators contribute to ^ — >■ 
7 + W at the tree level, 

a) \^^Alt''W^''W^,Ou, ^^Alt^^B^^'B^A. X^n,Au''^W^''W,,Ou, hbAu'" B^'B.^Ou, 
h) XqqA^'^^Ql^^D'^QlOu, XuuAu'"Unl,D^^UnOu, XddAu''" Drj^D>^DrOu, 

c) ^qqA^'^^Ql^^QlO'^Ou, Xuu^u''''URl,URd^Ou. XddAu''" Dnj^Dnd^^Ou, (18 

d) ^qqAIt'^Ql^.Q^O^, X^^^AI'-Ur^.UrO^, X^o^Al'-Dny.DnO^, KoAI'^ B.^d'^O^ 
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For class a) contributions, we obtain 



where Xgg takes the values Xww sin^ 9w , A;,;, cos^ 9w, X^w sin^ 9w and Xbb cos^ 9w for the four 
operators in class a) in order, respectively. 
For classes b) and c) contributions, we have 

Br{V^^ + U,X) _ f^^^ A,^mlE,X' 



Br{V fi+fi-) J ^87T^a{AlY'^{ml - 2mvE^f-'^u ' 



For class d), contributions from the first three operators are give by 
BriV^^+U,X') _ f^^ A,,{ml-mvE,)E,X'^ 



Br{V^fi+fi-) J ^2T^'^a{KliY^-^{ml-2mvE^f-^u 

In the above A = Xqq or Add and A = Xqq or Add for quarkonia composed of down and 
up type of quarks, respectively. Similarly for A and A'. 

The fourth operator X'^^qK]^'^^ B^i,d^O^ in class d) also contributes ioV 7 + W and the 
decay width can be obtained by replacing A' with X[Q{eQqCos9w) in eq.( l2Ti) . 

It is interesting to note that for contributions from classes a), b) and c), du needs to be 
larger than 1 in order to have a finite width for — > 7 + W, while for the contributions from 
class d), dy needs to be larger than 2 to have a finite width. In our numerical analysis, we 
will let du to be larger than 2 for this case. Also note that the distributions of E^ for class 
a), classes b) and c), and class d) are different. This can be used to distinguish different 
contributions if enough data are accumulated. 

There are experimental constrains on 7 + invisible decays of T and J /ip with Br{J/i/j 



7 + invisible) < 1.4 x IQ-^ |30|, and Br{T{lS) ^ 7 + invisible) < 1.5 x 10"^ |3l|, l32|. 
Combining the formula obtained above for unparticle contributions, one can set constraints 
on unparticle interactions. The results are shown in Fig. |H We find that the present upper 
bounds on T{J/ip) — 7 + invisible decay widths do not give strong constraints on the 
unparticle interactions. Improved bounds can provide more information. 

We comment that there are also contributions to Z +U from class a) operators. 
The decay width is given by 

dV _Ad^sm\29^) E^mz /^m| - 2mz^y" ^2 ^^2) 



dE^ 37r2 (m| - 2mzE^Y \ A, 



Y(1S) -> y + invisible 
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FIG. 4: Bounds on parameter space of A vs. Ku allowing unparticle decay mode to satm'ate the 
experimental data for 7 + invisible of T and J/V' decays. From top to bottom, the lines represent 
the cases for du = 1.3, 1.5, 2.1, 2.5 and for the bottom panel, du only takes 2.1 and 2.5 . 

where A can be any of \ww,bb, ^ww,bb- It would be interesting to see if strong constraints can 
be obtained for unparticle interactions when LEP data are analyzed for Z — > 7 + invisible. 
Summary 

If unparticles exist they must interact weakly with particles. The direct signature of 
unparticles will be in the form of missing energy in decays and collisions of particles. In this 
paper we have studied constraints on unparticle interactions using totally invisible decay 
modes of Z, vector quarkonia V and neutrinos. There are several operators which can 
contribute to these decays. 

Two operators with couplings and A^,^ contribute to Z ^U. Numerically the bound 
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on Im{X'i^f^) is about 5.5 times stronger than A'^q for given du and A^. The constraints are 
very sensitive to the dimension parameter du. If du is close to 1, with A'^^q = 1 and du = 1.3, 
Au needs to be larger than lO'* TeV, but Au can be as low as one TeV for du = 2. If by 
some means that the scale Ak is known, for example Ay = 1 TeV, we have the upper bounds 
Ko — 0.049 and 0.10 for the cases du=1.3 and 1.5, respectively. 

Several operators contribute to V ^ U, including the operator with coupling X[q and 
additional ones ^qq^j j^j^. There is experimental upper bound from T invisible. We find 
that the constraints on A^q are weaker than the constraint on Agg by a factor of 5.6. 
The constraints are again sensitive to du- If du is close to 1, with Agg = 1 and du = 1.3, 
Au needs to be larger than 2 x 10^ TeV, but Au can be as low as 400 GeV for du = 2. If 
the scale Au is set to be 1 TeV, we find the upper bounds Xqq dj:! = 0.025 and 0.082 for the 
cases du=l-3 and 1.5 respectively. The constraint on A'^q obtained from T — > W is weaker 
than that obtained from Z — W by a factor of 1.49(m^/mx)'^"~^, which is in the range 
2.9 ~ 14.4 when du takes value in the range of 1.3 ~ 2.0. 

There is one operator which can induce neutrino to U decay. Strong constraint could be 
obtained using constraint on r/m obtained from solar neutrino data. If one sets Au to be 
less than the Planck scale mp = 1.22 x 10^^ GeV, the dimension du must be bigger than 
1.6. 

We also studied the possibility of using V{Z) •y+U to constrain unparticle interactions. 
We find that present experimental upper bounds for T(J/^/') — * 7 + invisible does not give 
strong bounds on unparticle interactions. 
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